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Abstract 

We prove uniform Morrey-Campanato estimates for Helmholtz equations in the 
case of two unbounded inhomogeneous media separated by an interface. They imply 
weighted L 2 -estimates for the solution. We also prove a uniform L 2 -estimate without 
weight for the trace of the solution on the interface. 
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1 Introduction 

We consider the following Helmholtz equation 

ieu + Au + n(x)u = f(x), x G R d . (1.1) 

We assume that the refraction index n is nonnegative and discontinuous at 
the interface between two unbounded inhomogeneous media fi_|_, f2_ such that 
VL+ U TT = Tl~ U Q- = R d . We write 



n(x) 



n + (x) if x G S7_|_ 
n_(x) if x G 



We study here the limiting absorption principle, i.e. the limit when e — > + 
of equation (|l.l|l : our goal is to prove bounds on u that are uniform in e. A 
related question is the statement of the Sommerfeld radiation condition for the 
following equation 

Au + n(x)u = f{x), x G K d . (1.2) 
This is an open problem in the full generality of our assumptions, stated below. 



In this paper, we prove uniform Morrey-Campanato type estimates for this 
equation, using a multiplier method borrowed from [Hj. These bounds encode 
in the optimal way the decay: |u(x)| ~ l/\x\^~ at infinity of the solution u. 
They imply weighted L 2 -estimates for the solution u. We also prove a uniform 
L 2 -estimate without weight for the trace of the solution on the interface, which 
states that u carries essentially no energy on this set. 

In order to state precisely our assumptions and results, we need the following 
notations. 

First, we use the Morrey-Campanato norm, defined for u G Lf oc , 

Mb* = supi [ \u{x)\ 2 dx, (1.3) 

fl>0 -K Jb(R) 



where B(R) denotes the ball of radius R. We also use the following dual norm 



n/b=E( 2 ' +1 / i/wi 2 ^) 



1/2 



(1.4) 



where C(j) = {x e R d /2^ < \x\ < 2^' +1 }. 
The duality is given by the easy estimate 



/ fudx\ < \\.fU 



u 




(1.5) 



Note that both norms B and B* are homogeneous in space. 



We also denote the radial and tangential derivatives by 



dr \x\ ' 



d x 



d_ 



x d 
\x\ dr 



On the other hand, we assume that the interface between the two media 
T = dil + = <9£1_ is a smooth surface (Lipschitz is enough). Let dj be the 
euclidian surface measure on T and v(x) be the unit normal vector at x G T 
directed from to Wc denote, for x e T, the jump 



Throughout this paper we will write Vn instead of Vtz+1q + + Vn_lji_, the 
derivative of n outside the interface. Similarly, dd.n will denote the partial 
derivative of n with respect to the Xd variable outside the interface. 

The key assumptions we make on the interface T and the refraction index n 
arc the following. We comment on these assumptions later. 
(HI) There is a > such that the d-th component of v satisfies 



(H2) [n](x) has the same sign for all x e T; the following notation will be 
convenient: let a = — if [n] is non-negative and a = + if [n] is non-positive. 
(H3) n e L°°, n > 0. 
(H4) 



[n](x) = n + (x) — n_(x). 



Vd{x) > a for all x G V. 



(1.6) 




(1.7) 



(H5) 




(1.8) 



(H6) fr+fo < 1. 



We are now ready to state our main theorem. 
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Theorem 1.1 For dimensions d > 3, assume (H1)-(H6). Then, the solution 
to the Helmholtz equation satisfies the following estimates : 




)TA%+I l ^dx+sup± / \u\ 2 da R 
Jr" \x\ r>o R 2 Js(R) 



+ 



[\[n}\ \u\ 2 d 1 + [ \d d n\ \u\ 2 dx < C\\f\\%, (1.9) 



Jr JR d 



where C is a constant depending only on a, [3\ and 0%. 

We would like to stress several aspects of our analysis. First, the homogene- 
ity of the estimates and assumptions makes this theorem compatible with the 
high frequencies. The scaling invariance plays a fundamental role in the high 
frequency limit of Helmholtz equations (see Benamou et al [5] , Castella et al |3] 
for the case of a regular index of refraction). 

We would like to point out several terms in the left-hand side. The last two 
say that, in principle, the energy is not trapped at the interface and it mainly 
radiates in the directions where ddn vanishes. On the other hand, the first two 
terms essentially assert that u and Vu belong to the optimal space B* , provided 
/ 6 B. This is the same Morrey-Campanato estimates as in the regular case. 

Let us now comment our assumptions on n. First, they allow some growth 
at infinity : n does not go to a constant at infinity. The hypotheses (H2), (H4) 
and (H5) can be understood as conditions on the trajectories of the geometrical 
optics. The condition (H4) implies the dispersion of these trajectories. The 
conditions (H2) and (H5) ensure that the energy goes from one side of the in- 
terface to the other. The condition (H5) involves both the interface and the 
index: it becomes a weaker assumption on the index when the interface is close 
to a hyperplane (a ~ 1). This type of assumptions is natural in the study of the 
high frequency limit where the link with Liouville's equations can be understood 
through Wigner transform (see L. Miller [5] for a refraction result in the case of 
a sharp interface for Schrodinger equation, E. Fouassier [H] for high frequency 
limit of Helmholtz equations with interface; for an account on high frequency 
limit for wave equations, see P.-L. Lions, T.Paul [Jj, and 0, for Helmholtz 
equations without interface). 

Taking [n] = in Theorem ll.il i.e. the case without interface, our results 
boil down to the uniform estimate proved by B. Perthame, L. Vega jU] when 
the refraction index satisfies assumptions (H3)-(H4) with /3x < 1. Under these 
assumptions, they also proved in |l()j . jll| an energy estimate saying that the 
energy \u\ 2 mainly radiates in the directions of the critical points of (where 
Hoc is given by n(x) — > n oc (-^), \x\ — > oo). In our case, the corresponding 
energy estimate corresponds to the last two terms. 

In the case of two unbounded media, similar results, but not scaling invariant, 
were obtained in previous papers. Eidus [2] first proved weighted-L 2 -estimates, 
and the L 2 -estimate on the trace of the solution on the interface, for a piecewise 
constant index of refraction. To do so, he assumed that the interface satisfied 
an extra "cone-like" shape condition, \x ■ v\ < C for x G T. Under the same 
assumption on the interface, Bo Zhang also proved inhomogeneous B*-B 
estimates when n is a long-range perturbation of a piecewise constant function. 
S. DeBievre, D.W. Pravica 0] proved weighted-L 2 -estimates in a very general 
context using Mourre's commutator method. In particular they considered the 
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case of an index that is smooth outside of a compact set in the x^-direction and 
only bounded in this compact set. 

Our proof is based on a multiplier method. Following B. Pcrthame, L. 
Vega 0, we use a combination of a Morawetz-type multiplier and of an ellip- 
tic multiplier. Following Eidus 0, we combine it with a multiplier specific to 
the case with an interface for which one direction plays a particular role (the 
x^-direction here). The first two multipliers allow us to control both Vm and 
u locally in L 2 by \\f\\g and we estimate the integral over the interface of \u\ 2 
using the third multiplier. 

The paper is organized as follows. In Sections 2 and 3, we present the 
basic multipliers and the particular choice we make here to prove Theorem ll.il 
Then, in Section 4, we give another estimate containing the trace of Vw on the 
interface. 



2 Basic identities 

Lemma 2.1 The solution to the Helmholtz equation satisfies the following 
four identities, for smooth real valued test functions tp, ip, 

ip(x)\Vu(x)\ 2 + — / Aip(x)\u(x)\ 2 + / ip(x)n(x)\u(x)\ 2 

2 Imd lad 



= TZe / f(x) V {x)u(x), (2.1) 



e I ip(x)\u(x)\ 2 — Zm / V<p(x) ■ \/u(x)u(x) = Im I f (x)(p(x)u(x) , (2.2) 



Vu(x) / A 2 tfj(x)\u(x)\ 2 + — / Vn(x)-\7if)(x)\u(x)\ 2 

4 jRd 2 J R d 

-elm I W(a;) • Vu(x)u(x) + - / \n]v{x) ■ W^(x)\u(x)\ 2 d-y(x) 

= -TZe I f(x)(Vip(x)-Vu(x) + ^Aip(x)u(x)), (2.3) 



\ [ [n]i' d (x)\u(x)\ 2 d-f(x) + i I d d n(x)\u(x)\ 2 

Z JV Z JR d 

= —Tie / f(x)ddu(x) — elm I u{x)ddu(x). (2.4) 

JR d JR d 

Proof: The identities l|2.1|l and l|2.2ll are obtained by multiplying the Helmholtz 
equation Hl.lfl by ipu and then taking the real and imaginary parts. The iden- 
tity H2.3|) is obtained using the Morawetz-type multiplier Vijj(x) ■ Vu(x) + 
7}Aip(x)u(x) and taking the real part. To get the last identity H2.4J) . we use 
the multiplier ddU and take the real part. □ 
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Lemma 2.2 The solution to the Helmholtz equation J^l.l}) satisfies the following 
estimate 

\[n]\\u{x)\ 2 d 1 {x) + - f (d d n) a \u\ 2 
r a JR d 

<2 [ \fd d u\ + 2e\lm [ ud d u\ + - [ (d d n)_ a \u\ 2 . (2.5) 
jR d Jm d a JR d 

Proof: This lemma follows directly from the identity (|2.4|l using the hypothesis 
(HI). □ 

3 Proof of theorem 11.11 

The following proof in which all the details are included for the convenience 
of the reader is essentially adapted from j^J, apart from the treatment of the 
terms involving the interface and the partial derivative with respect to the x d - 
direction, which is the key difficulty of this paper. 

We derive the proof of theorem II . II from the above identities. We make the 
following choice of test functions ip and (p, for R > 0, 



Vip{x) 



x/R for \x\ < R 
x/\x\ for \x\ > R, 



1/2R for |x| < R 
for \x\ > R. 



tp(x) 

We also need the following calculations (in the distributional sense) 



D^(x) = \ 5ij/R tor\x\<R 
y 1 (5ij\x\ 2 — XiXj)/\x\ 3 for \x\ > R, 



d/R for |x| < R 

for \x\ > R, 



and the inequality 

1 f d 1 f 

- / vk{2 V - A^) > — \ vda R for v > 0. (3.1) 

4.V 4-R Js(R) 

We add the identity (|2.1J) to (|2.3(l . which gives, using the inequality l|3.1J) . 
for the previous choice of "0, tp, 



2R 



'B(R) 



2R 



— |Vu(x)r + — / n\u{x)\ 2 + n |Vo(a) 



B(R) 

d-l 
4i? 2 



\x\>R 



x ^ 



S{R) 



\u\ 2 da R + - / (ViP(x)-Vn(x)) + \u(x)f 



K*)| 



<C \f(x)\^Jl+C \f(x)\\Vu(x)\ + - (V^(x)-Vn(x))_Kx)| 



+ Ce I |u(a)||V«(a?)H-^ / |N| \v ■ ViP\\u(x)\ 2 d 1 . 
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Then, we use the inequality 1)2. 5[) . together with the bound |V^| < 1 to 
estimate the trace term 



1 / |Vu(a;)| 2 + i / | U (x)| 2 + i / (x • Vn(s)) + K*)| 2 (3.2) 

^B(fl) n JB(R) n JB(R) 

f |V T u(a;)| 2 d- 1 f ., , 

•/|£c|>,R Fl ZJX JS(R) 

<C I \f(x)\ l -^ + C [ \f(x)\\Vu(x)\+Ce f |ti(s)||V«(aO| (3.4) 

JR d Fl JR d JR d 

+ [ (V4>(x)-Vn(x))-\u(x)\ 2 + - [ (d d n{x)) a \u(x)\ 2 . 

JR d a JR d 

(3.5) 

Our task is to estimate the terms in the right-hand side of the last inequal- 
ity. We separate them in three types : those containing the source /, those 
containing e and those containing the index n. 

We begin by the two terms containing /. Using the duality estimate (|1.5|) . we 
get for all S > 

/ |/||Vtt| < Iivub.11/H5 

JR d 

< 5\\Vu\\%+Cs\\f\\l. 

For the second term, we use a Cauchy-Schwarz inequality to obtain, for all 
6 > 



< 



CO') Ft/ \ Jc(j) 



1/2 / ^ \ 1/2 

' 2 



supij/ H 2 ^) W2' / l/l 
< <Ssup4 / M 2 da fl + CW(/) 2 



?>0 ^ Js(fl) 

Next, we turn to the term containing e. First, using the identities l|2.2|) and 
(|2.1() with tp = 1, one can notice 



M 2 < / IM 

JR d 

s I \Vu\ 2 < \\n\\ L ~ f \fu\. 

JR d jR d 

Hence, using again the duality estimate l|1.5|) . we obtain for all 5 > 
e I \u\\Vu\ < C\\u\\s4f\\ B 

JR d 

< S\\Vu\\% m + C s \\f\\%. 
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The terms with n remain. We write 

,2 ( x ■ Vn(a;))- 



/ (Vt/;- Vn)_| U | 2 < V / 



< 



Rtl 



i?>0 « Js(fl) / ^ C(j) 



Hence. 



Similarly, we get 



(Wip(x) • Vn(x))_Ka;)| 2 < /?i||u[ 



- / (d d n(*))<rK*)l" <f3 2 \\uf g „ 



a 

Putting all these estimates together gives for all S > 

_L /" | Vw (x)| 2 + I / | U (x)| 2 + I / (rr ■ Vn(x)) + | U (x)| 2 
k Jb(r) k Jb(r) u Jb(r) 

V T u(x)\ 2 , d-1 ( 



|x|>R 1*1 2 ^ 2 JS(fl) 



< (/3 1 + /3 2 + <5)H||,+C 5 ||/|||. 

Hence, choosing 5 small enough (depending on /3i + ,$2) and taking the 
supremum with respect to R, we obtain 

UVulljL + \\n^u\\%, + I ^ + ||(* • Vn)f u\\% (3.6) 



K d 1*1 



+ sup-^/ M 2 ^, < C||/|||. (3.7) 

To end up the proof, it remains to estimate the two last terms in the left-hand 
side of l|1.9|l . We use the inequality Q2.5JI and the previous bounds 



(d d n) a \u\ 



[n]\ \u(x)\ d'y + (d d n)- a \u\ 
<C I \fd d u\+Ce j \u\Vu\ + C- 

<c(||/||| + H|.) 

Using l|3.6|l we obtain 

n]| \u{x)\ 2 d 1 + I (d d n)^ a \u\ 2 <C\\f\\l 



r 



This ends the proof. 

4 Another trace estimate 

In this section, we assume that the interface T is a hyperplane T = {x d = 0}. 
Then we prove the following extra uniform estimate 
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Theorem 4.1 Under the assumptions of Theorem \l.l\ if we assume moreover 
that r = {xd = 0} and there exists (3 > such that (x) 1+ ^\ V^ml <E L°° , then 

Jj[n]\ \Vu(x)\ 2 dx' <c(\\f\\% + \\V x ,f\\%) (4.1) 

where C is a constant depending only on a, (3% and IKa;) 1 V x 'n||L~. 

To prove this theorem, we will need the following identity 
Lemma 4.1 The solution to the Helmholtz equation satisfies 

\ f [n] \V x ,u(x)\ 2 dx' -\f[n] \d d u{x)\ 2 dx' - \ f [n 2 ] \u(x)\ 2 dx' 
•J r j r •/ r 

dd{n 2 {x))\u\ 2 — elm I n{x)u{x)ddu{x) 

JR d 

= TZe [ n(x)f(x)d d u(x). (4.2) 

Proof: This identity is obtained by multiplying the equation 11.1|) by nddU and 
taking the real part. □ 
Proof of Theorem 14.11 

Firstly, note that V^n u E B: 



x) 2 - 2 P\u\ 2 ) 1/2 



\\V x ,nu\\ B < ]T(2^ / ( 

E( 2J+ 7 < 



< 

j<0 



x) 2 - 2 ^\u\ 2 ) 1/2 



j>0 

but, using H1.9J) . we have, for S > 0, 

\\(x)-^ 5 u\\ L2 <C\\u\\^<C\\f\\ 6 , 

hence 

\\v x ,nu\\ B < c\\f\\ B j2^+cj2^ J " /2 u x y L2 ^ u \\^ 

j<0 j>0 

< cil/IU- 

So, we can apply the estimate (|1.9|l to V x >u. Indeed, V X 'U satisfies the following 
Helmholtz equation 

ieW x 'U + AV X 'U + n{x)V X 'U = V x >f(x) + V x /n u, (4-3) 

so V T 'ti satisfies 



\[n]\ \VM 2 dx' <c(\\(V x ,ff B + \\f\\ B ). 
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Then, we use (14 .2[) to get 

/ |[»]| \d d u\ 2 dx'< f \[n}\ \V x ,u\ 2 dx' + f \[n 2 }\ \u\ 2 dx' + f \d d (n 2 )\ \u\ 2 
Jr Jr Jr J~R d 

+ e f n\ud d u\ + f n\fd d u\. (4.4) 

jR d JR d 

Since n G L°° , we have 

J \[n 2 ]\ \u\ 2 dx' < 2||n|| oo y < C||/||| 

and 

/ |0«i(n 2 )|M a < 2||n||oo / |9 d n|H 2 <(7||/|||. 

JR d JS. d 

Hence, estimating the two last terms in the rigth-hand side of (|4.4|l as before 
and using we get 

J \[n}\ |^| 2 dx'<c(||/||| + ||V^/||s), 

and thus the theorem is proved. 
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